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Abstract 

The dynamical equation of quantum mechanics are rewritten in form of 
dynamical equations for the measurable, positive marginal distribution of the 
shifted, rotated and squeezed quadrature introduced in the so called "sym- 
plectic tomography". Then the possibility of a purely classical description of 
a quantum system as well as a reinterpretation of the quantum measurement 
theory is discussed and a comparision with the well known quasi-probabilities 
approach is given. Furthermore, an analysis of the properties of this marginal 
distribution, which contains all the quantum information, is performed in the 
framework of classical probability theory. Finally examples of harmonic oscil- 
lator's states dynamics are treated. 
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1 Introduction 



" Schrodinger made no secret of his intention to substitute simple classical pictures 

for the strange conceptions of quantum mechanics, for whose abstract character he 
expressed deep aversion" . It is clear from this commentary of Rosenfeld H] that from 
the early days of quantum theory there has been a permanent wish to understand 
quantum mechanics in terms of classical probabilities. However, due to the Heisem- 
berg [[J and Schrodinger- Robertson ||, [|J uncertainty relation for the position and 
momentum in quantum systems, does not exist joint distribution function in the 
phase space. This leads to the introduction of the so called quasi-probability distri- 
butions, such as Wigner function ||, Husimi Q-function || and Glauber-Sudarshan 
P-function [[5J, |§; later on unified into one-parametric family |§. Furthermore, in 
order to get a bridge between quantum and classical physics, Madelung fl(| already 



observed that the modulus and the phase of the wave function obey the hydrodinam- 
ical classical equations, and along this line the stochastic quantization scheme has 



been suggested by Nelson [ I I] to link the classical stochastic mechanics formalism 



with the quantum mechanical basic entities, such as wave function and propagator. 



In some sense, also the hidden variables [12] was proposed to relate the quantum 



processes to the classical ones. Nevertheless, up to date does not exsist a formalism 
which consistently connects the "two worlds". 

The discussed quasi-probabilities illuminated the similarities and the differences 
between classical and quantum considerations, and they are widely used as instru- 
ments for calculations in quantum theory [T^|, [T^|. However, they cannot play the 
role of classical distributions since, for example, the Wigner function and the P- 
function may have negative values. Althought the Q-function is always positive and 
normalized, it does not describe measurable distributions of concrete physical vari- 
ables. 



Recentely, after J. Bertrand and P. Bertrand |k| made the first attempt to apply 
the tomographic principle to phase space distributions, Vogel and Risken [16|, using 
the formalism of Ref. 0, established an integral relation between the Wigner function 
and the marginal distribution for the measurable homodyne output variable which 
represents a rotated quadrature of the electromagnetic field. This result gives the 
possibility of "measuring" the quantum state, and it is referred as optical homodyne 
tomography |L7 . 



In Ref. [|18j a symplectic tomography procedure was suggested to obtain the 
Wigner function by measuring the marginal distribution for a shifted, rotated and 



squeezed quadrature, which depends on extra parameters. In Ref. the formalism 
of Ref. |16[] was formulated in invariant form, relating the homodyne output distribu- 
tion directly to the density operator. Then, in Ref. |2(| the symplectic tomography 
formalism was also formulated in this invariant form and it was extended to the mul- 
timode case. Thus, due to the introduction of quantum tomography procedure the 
real positive marginal distribution for measurable observables, such as rotated shifted 
and squeezed quadratures, turned out to determine completely the quantum states. 

The aim of the present work is to formulate the standard quantum dynamics in 
terms of the classical marginal distribution of the measurable shifted, rotated and 
squeezed quadrature components, used in the symplectic tomography scheme. Thus 
we obtain an alternative formulation of the quantum system evolution in terms of 
evolution of a real and positive distribution function for measurable physical observ- 
ables. We will show the connection of such a "classical" probability evolution with 
the evolution of the above discussed quasi-probability distributions. Preliminarly, the 
approach was shortly presented in Ref. |2l[ . 

Examples relative to states of harmonic oscillator and free motion will be consid- 
ered in the frame of the given formulation of quantum mechanics as well as oscillator 
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with friction and driven terms included. 



2 Density operator and distribution for shifted ro- 
tated and squeezed quadrature 



In Ref. [18] it was introduced an operator X as the generic linear combination of the 



position q and momentum p (% = 1) 

X = fiq + up + 5, (1) 

which depends upon three real parameters /x, u, 6 and, due to its hermiticity, is 
a measurable observable. Thus, the marginal distribution, defined as the Fourier 
transform of the characteristic function 



w 



(JC./x, !/,*) = / dke~ ikx (e ikx ), (2) 



depends itself upon the parameters /j, u, 5, and it is normalized with respect to the 
X variable 

J dXw{X,n,v,5) = 1. (3) 

Furthermore, it was shown |18| that this marginal distribution is related to the state 
of the quantum system, expressed in terms of its Wigner function W(q,p), as follows 

<X,^S) = I e-^ x ^-^W( q ,p)^^. (4) 

Eq. @ shows that w is a function of the difference X — 5 = x, so that it can be 
rewritten as 

w(x,n,v) = I e- lk ^-^W(q,p)^P$. (5) 



(2vr) 2 

This formula can be inverted and the Wigner function of the state can be expressed 



in terms of the marginal distribution [18 



W(q,p) = (2tt) z w f (z, -zq, -zp), (6) 



where Wp(z, a, b) is the Fourier component of the marginal distribution (||) taken with 
respect to the variables x, //, v, i.e. 

w F (z, a, b) = — [ w(x, /i, v)e- l{xz+ ^ a+vb) dxd^dv. (7) 
{2ix) A J 

Hence, it was shown that the quantum state could be described by the positive 
classical marginal distribution for the squeezed, rotated and shifted quadrature which 
could be considered as a classical probability associated to a stochastic variable x and 
depending also on parameters. 

In the case of only rotated quadrature, p = cos0 and v = sin0, the usual optical 



tomography formula of Ref . JT6[ , gives the same possibility through the Radon trans- 
form instead of the Fourier transform. This is, in fact, a partial case of the symplectic 
transformation of quadrature since the rotation group is a subgroup of the symplectic 
group ISp(2,R) whose parameters are used to describe the transformation (|l]). 

In Ref. pO] an invariant form connecting directly the marginal distribution 



w(x, p, v) and the density operator was found 

p = J dxdfidv w(x, /i, v)K^ v , (8) 
where the kernel operator has the form 

The formulae (^) and ([8p of symplectic tomography show that there exist an invertible 
map between the quantum states described by the set of nonnegative and normalized 
hermitian density operators p and the set of positive, normalized marginal distribu- 
tions (" classical" ones) for the measurable shifted, rotated and squeezed quadratures. 
So, the information contained in the marginal distribution is the same which is con- 
tained in the density operator; and due to this, one could represent the quantum 
dynamics in terms of evolution of the marginal probability. Really, the fact that 



Kii,v depends also on the z variable (i.e. each Fourier component gives a selfconsis- 
tent kernel) shows the overcompleteness of information achievable by measuring the 
observable of Eq. (fj) . 

The definition of the marginal distribution function w(x,fi,u) might be alterna- 
tively given in terms of the eigenstates of the operator x = X — 5 

x\x) = x\x) (10) 
which can be obtained from the position eigenstates 

q\q)=q\q) (11) 
by the action of the unitary operator S 

\x) = S\q) (12) 

which represents the composition of simple operations such as rotation and squeezing, 
i.e. it satisfies the requierement 

S^qS = jiq + isp. (13) 



It is worth to remark, about this transformation, that there exist a costraint due 
to the commutation relation between the observable (|I|) and its canonical conjugate, 
i.e. if one introduce the observable 

p = ^ + u'p + S' (14) 

the matrix 



must satisfy the relation 



Kak T = a- 0=[\ J). (16) 



Then, the marginal distribution is the diagonal matrix element of the density 
operator in the transformed basis fllOf) 



w(x, p, v) = (x\p\x) = Tr{p\x)(x\} (17) 

or it is the diagonal matrix element in position representation of the transformed 
density operator 

w(x,p,u) = (q\&pS\q) = Tr{& pS\q)(q\}. (18) 



The form of the shifted and squeezed operator S is well known |23| . Choosing the 
parameters p = cos</> and v = sin0, the operator S gives the marginal distribution 
for the homodyne output of Ref . |Tj| . In the case of p = 1 and v — the marginal 



distribution is that for quadrature q, i.e. w(q, 1,0) = p(q,q) = (q\p\q), while in the 
case of p = and v = 1 the marginal distribution is that for the other quadrature p, 
i.e. w(p, 0, 1) = p(p,p) = (p\p\p). 

3 Quantum evolution as a classical process 

We now derive the evolution equation for the marginal distribution function w using 
the invariant form of the connection between the marginal distribution and the density 
operator given by the formula (§). Then, from the equation of motion for the density 
operator which includes the interaction with environment x{p) 

d t p=-i[H,p]+ X (p), (19) 
we obtain the evolution equation for the marginal distribution in the form 

dxdpdv {w(x, p, v, t)K^ u + w(x, p, u, t)/ M ^| = x (J dxdpdv w(x, p, u, t)K^ u 

(20) 

in which the known Hamiltonian determines the kernel through the commutator 

I^ = i[H,k^\, (21) 



while the r.h.s. is functionally dependent on the marginal distribution. The obtained 
integral-operator equation can be reduced to an integro-differential equation for the 
function w in some cases. Let us consider at first the situation in which x(p) = 0, 
the opposite situation will be discussed later. Then, we represent the kernel operator 
7 M)iy in normal order form (i.e. all the momentum operators on the left side and the 
position ones on the right side) containing the operator K^ v as follow 

: 4„ := K(p) : : V(q) (22) 

where TZ(p) and V(q) are, finite or infinite operator polynomials (depending also 
on the parameters /i and v) determined by the Hamiltonian. Then we calculate the 
matrix elements of the operator equation (Effi ) between the states (p\ and \ q) obtaining 

J dxdfidv {w(x, /i, v, t) + w(x, /i, v, t)H(p)V(q)} (p\ : K^ u : \q) = . (23) 

If we suppose to write 

K(p)V(q) = n(p, ?) = EE c n,m(z, V, v)V n q m , (24) 

n m 

due to the particular form of the kernel in Eq. (^), the Eq. (££3) can be rewritten as 
J dxd[idv \w{x, fx, v, t) + w(x, /i, u, t)Tt(p, g)| (p\ : K^ u : \q) — , (25) 
where p, q are operators of the form 

while in the space of variables x, /x, z/ should be intended as the derivative with 
respect to x, i.e. 

2 " 4 (27) 

and when it appears in the denominator is understood as an integral operator. Fur- 



thermore the right arrow over II means that, with respect to the order of Eq. (|24) , 



the operators p and q act on the right, i.e. on (p\ : K^ v : \q). Under the hypothesis of 
regularity of w on the boundaries, we can perform integrations by parts in Eq. (|25f) 
disregarding the surface terms, to get 

J dxdfidv {w(x, fj,, v, t) + w(x, /i, u, t) tl (p, g)| (p\ : K^ u : \q) — , (28) 

where now LT means that the operators p, q 

'-Hl-s-)- '-Hi-?)' < 29 > 

act on the left, i.e. on the product of coefficients c n ^ m {— z, fi, v) with the marginal 
distribution w. Finally, using the completness property of the Fourier exponents 
given by (p\ : K^ v : |g) we arrive at the following equation of motion for the marginal 
distribution function 

d t w + wtT (p, q) = . (30) 

Let us consider the important example of the motion of the particle in a potential 
with the Hamiltonian 

H = ^ + V(q); (31) 



then the described procedure of calculating the normal order kernel ( P2|) gives the 
following form of the quantum dynamics in terms of a Fokker-Planck-like equation 
for the marginal distribution 



d 

w — Utt w ~ 1 
ov 



\ T (-Id .v d \ T / -1 d .v d 



w = (32) 



d/dxdfj, 2 dx J \d/dxdpL 2dx / 

which in the general case is an integro-differential equation. It is worth to remark that 
considering the quadrature X of Eq. (fj) to be dimensionless, the Planck constant H, 
should appears in the Eq. ([32]) to multiply the fisrt two terms. As a consequence it 
is clear that the equation, even if classical-like, gives a quantum description of the 
system evolution (as the Schrodinger equation). 

9 



Thus given a Hamiltonian of the form fl3~T|) we can study the quantum evolution of 
the system writing down a Fokker-Planck-like equation for the marginal distribution. 
Solving this one for a given initial positive and normalized marginal distribution we 
can obtain the quantum density operator p(t) according to Eq. (H). Conceptually it 
means that we can discuss the system quantum evolution considering classical, real 
positive and normalized distributions for the measurable variable X which is shifted, 
rotated and squeezed quadrature. The distribution function which depends on extra 
parameters obeys a classical equation which preserves the normalization condition 
of the distribution. In this sense we always can reduce the quantum behaviour of 
the system to the classical behaviour of the marginal distribution. Of course, this 
statement respects the uncertainty relation because the measurable marginal distri- 
bution is the distribution for one observable. That is the essential difference (de- 
spite of some similarities) of the introduced marginal distribution from the discussed 
quasi-distributions, including the real positive Q-function, which depend on the two 
variables of the phase space and are normalized with respect to these variables. We 
would point out that we do not derive quantum mechanics from classical stochastic 
mechanics, i.e. we do not quantize any classical stochastic process, our result is to 
present the quantum dynamics equations as classical ones, and in doing this we need 
not only the classical Hamiltonian but also its quantum counterpart. 

4 Examples 

Let us choose as system to study a driven harmonic oscillator of unit mass with an 
hamiltonian of the type 

H = ^+^-fq, (33) 
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then from Eq. (|32|) immediately follows 



d d d 

w — u—w + lu 2 v—w + fv—w = . (34) 
av a/x ox 



Below we cosider solutions of some special cases of Eq. ([54]) , while the solution for 
the complete equation will be given in a next section by using a propagator method. 

4.1 Free Motion 

For the free motion, u — f — 0, the evolution equation (|3~4|) becomes the first order 
partial differential equation 

w — Utt w = ) (35) 
av 

and it has a gaussian solution of the form 

„(*,„,,«)= * expj-^} ,36) 



where the dispersion of the observable x depends on time and parameters as follow 

*x(t) = ^[/i 2 (l + t 2 ) +is> + 2fjust]. (37) 
The initial condition corresponds to the marginal distribution of the ground state of 



an artificial harmonic oscillator calculated from the respective Wigner function [0 . 
4.2 Harmonic Oscillator 

For the simple harmonic oscillator with frequency u — 1, we have / = then Eq. 
([34] ) becomes 



d d 

w — a— w + v— w = . (38) 
av a ijl 

If we consider the first excited state of the harmonic oscillator, we know the 
Wigner function 



Wi(q,p) = -2(1 - 2q 2 - 2p 2 ) exp[-q 2 - p 2 ]. (39) 
11 



It results time independent due to the stationarity of the state, but for small g and 
p it becomes negative while the solution of Eq. ( j38f) 

2 _s f x 2 ) 

w x (x,ii,v,t) = —j=[n 2 + v 2 \ 2 a; 2 exp <^ — — -\ (40) 

V 71 " { V z + is z ) 

is itself time independent, but everywhere positive. 

Indeed, a time evolution is present explicitly in the coherent state, whose Wigner 
function is given by 

W c (q,p) = 2exp{-g 2 - gg - p 2 -pi + 2(qq + pp ) cost - (pq - qp ) sint} (41) 

where go an d Po are the initial values of position and momentum. For the same state 
the marginal distribution shows a more complicate evolution 

w c (x,fi,v,t) = -^[fi 2 + u 2 ]-^ (42) 



7T 

( 1 

x exp < -ql -pi - — + 2-(p cost - g sint) 



x exp 



p? + v 2 



a r 
—x + go(/isint + 2/ cost) +po( z/ sint — /icost) 



It is also interesting to consider the comparison between the Wigner function and 
the marginal probability for non-classical states of the harmonic oscillator, such as 
female cat state defined as 



|a_) = iV_(|a) - I - a)), a = 2- 1 / 2 (g + ip ) (43) 

with 

f exp[(gg+p 2 )/2] ^ 
- l4sinh[(g 2 +p 2 )/2]J 1 j 

and for which the Wigner function assumes the following form 

W-(q,p) = 2N 2 e~ q2 ~ p2 { e~ q2 ~ p ° cosh[2(gg + pp ) cost + 2(gp - pqo) sint] 

- cos[2(gp -pqo) cost - 2(gg + pp ) sint]}. (45) 
12 



The corresponding marginal distribution is 



with 



and 



w_(x, [Aji/jt) = N^_[wa(x, fj,,u,t) — wb(x, fl,U,t) 

- w* B (x,fj,, v,t) + w A (-x,fx, u,t)] (46) 



w A (x, n, u, i) 



71 



X 



X 



x exp^-g -p + 2-(p cost - q smt) 



(47) 



X eX P 1 2 i 2 



—x + q (/j,smt + vcost) + p (jysmt — ficost) 



w B {x,fi,u,t) = —=[fi 2 + u 



7T 



(4J 



I 0l> ijC 

x exp < — 2i— (q cost + p sint) 



x exp 



■1 



-i— x + qo{ficost — v sin t) + po(fismt + ucost) 



^/j? + v 2 

The presented examples show that, for the evolution of the state of a quantum 
system, one could always associate the evolution of the probability density for the 
random classical variable X which obeys "classical" Fokker-Planck-like equation, and 
this probability density contains the same information (about a quantum system) 
which is contained in any quasi-distribution function. But the probability density has 
the advantage to behave completly as the usual classical one. The physical meaning 
of the "classical" random variable X is transparent, it is considered as the position in 
an ensemble of shifted, rotated and scaled rest frames in the classical phase space of 
the system under study. We could remark that for non normalized quantum states, 
like the states with fixed momentum (De Broglie wave) or with fixed position, the 
introduced map in Eq. (ffl) may be preserved. In this context the plane wave states of 
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free motion have the marginal distribution corresponding to the classical white noise 
as we shall see below. 



4.3 Squeezed Coherent States 

Here we will consider the marginal distribution w(x,fi, v) for the squeezed coherent 
states of the harmonic oscillator. Since the Wigner function of these pure gaussian 
states may be represented in the form pi 



W a (q,p) = 2exp 
where m is the dispersion matrix 



-\{ P -p{tu-m m-M p q _*® 



m 



a p (t) a pq {t) 

<Tpq(t) <Tq{t) 



and p(t) and q(t) are the mean values of the quadratures 



q(t) = y/m{ 



ae 



p(t) = >/23f( 



ae 



(49) 



(50) 



(51) 



The variances in eq. (pO] ) are given by 

1 



a p (t) 
a p (t) 



s cos 2 1 H — sin 2 1 



cos 2 1 + s sin 2 1 



2 

2 Vs 

-Is ) sin t cos t , 

2 V s; 



w(x, /i, V, t) 



1 f Jx - fiq(t) - vp(t)} 2 

exp < 



2ira x (t) 



2a x (t) 



where 



a x{t) = -[(i 2 a q {t) + v 2 a p (t) + 2/iva pq (t)] . 



(52) 
(53) 
(54) 



with s the squeezing parameter. Using Eq. (|^) in the formulae (|6]) and (0), we 
obtain for the marginal distribution the expression 



(55) 



(56) 
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Let us now take the limit s — > 0, this means that our marginal distribution becomes 
a delta function 



lim w = 5 (x — (fiq Q + up ) cost — (fip Q — uq ) sint) , (57) 

and as a consequence its spectrum will be constant and equal to unity for each values 
of the variable conjugate to X, thus it will correspond to the white noise spectrum. 
On the other hand, the nonnormalized quantum states, like the states with fixed 
momentum (De Broglie wave) or with fixed position, have a marginal distribution 
normalized and everywhere equal to one. Thus plane wave states of free motion 
correspond to the classical white noise distribution. 

5 Evolution in the Presence of Environmental In- 
teraction 

When a system is coupled with the "rest of Universe" the time evolution of the 
density operator is no longer unitary, and to treat the problem at quantum level, 
one needs of some approximations; usually the starting point is a simple system as 
an harmonic oscillator which linearly interacts with a bath idealized as an infinity 
of other harmonic oscillators, then the (master) equation for the density operator 



becomes [27 



P = ~i[a ] a, p] + x(p) 

1 t + + 

x(p) — —(rt + l)(2apa'—a'ap — pa'a) 

+ —n(2a^pa — aa) p — paa^) (58) 

where 7 is the damping constant characterizing the relaxation time of the system, 
a, are the boson operators of the system and n is the number of the thermal 
excitation of the bath. Using Eq. ([H|) in the interaction picture and performing step 
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by step the same procedure that leads to Eq. (|30D, one may describes the damped 
evolution by means of 



w 



7 



9 d l.o 2 . d 2 



w 



(59) 



where we have assumed for simplicity n = 0, a situation common in quantum optical 
systems. In Eq. ( |59"|) we recognize the Fokker-Planck equation where the diffusion 
term is given by the proper stochastic term while the drift by the parameters (the 
factor 2 can be eliminated by a simple transformation w = we 7 '). The solution of 
Eq. (|59"D, with coherent initial excitation q , p , is 



exp 



-7*/2 



/J 2 + V 1 



(60) 



which is exactely the Fourier transform of the Wigner function for the damped har- 
monic oscillator given by |24 



W(q,p) = 2 exp -(q - g e' 7t/2 ) 2 - (p - R,e' 7t/2 ) 



- 7 t/2N2' 



(61) 



This is a proof that the developed formalism is consistent also in the case of open 
quantum systems. 

6 Quantum measurements and classical measure- 
ments 

In this section we will discuss the concept of quantum measurements in the frame 
of the developed approach. It is a well known steatment |28|], |29| that quantum 
mechanics suffers from an inconsistence in the sense that it needs, for its understand- 
ing, of a classical device measuring quantum observables. Due to this the theory of 
measurements suppose that there exist two worlds: the classical one and the quan- 
tum one. Of course in the classical world the measurements of classical observables 
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are produced by classical devices. In quantum world the measurements of quantum 
observables are produced by classical devices too. Due to this the theory of quantum 
measurements is considered as something very specifically different from the classical 
measurements. 



Recentely it has been proposed some schemes IpOR , [3jJ to resolve the dichotomy 
between the measured microsystem and the measuring macroapparatus, however it is 
phsycologically accepted that to understand the physical meaning of a measurement 
in classical world is much easier than to understand the physical meaning of an 
analogous measurement in quantum world. 

Our aim is to show that in fact all the roots of difficulties of quantum measure- 
ments are present in the classical measurements as well. Using the invertible map, 
of the quantum states (both normalized and nonnormalized) and classical states (de- 
scribed by classical distributions-generalized functions), given by Eq. (f|) we could 
conclude that the complete information about a quantum state is obtained from 
purely classical measurements of the position of a particle, made by classical devices 
in each reference frame of the ensamble of the classical reference frames, which are 
shifted, scaled and rotated in the classical phase space. 

These measurements do not need of any quantum language, if we know how to 
produce, in the classical world (using the notion of classical position and momentum), 
reference frames in the classical phase space differing from each other by rotation, 
scaling and shifting of the axis of the reference frame and how to measure only the 
position of the particle from the point of view of these different reference frames. 
Thus, knowing how to obtain the classical marginal distribution function w(x,fi,i>) 
which depends on the parameters /z, u, 6, labeling each reference frame in the classical 
phase space, we reconstruct through the map (H) the quantum density operator. 

By this approach, we avoid the unpleasant paradox of quantum world which 



17 



needs for its explanation measurements by a classical apparatus. Nevertheless all the 
difficulties of the quantum approach continue to be present, but in a different classical 
form. In fact, if we consider for example the notion of wave function collapse p2 |, it 



is displaced in the classical framework, since if we idealize the measuring apparatus 
as a bath with which the system interacts |p3|| , then a reduction of the probability 



distribution (as our marginal distribution) occours as soon as we "pick" a value (hence 
a trajectory) of the classical stochastic process associated to the observable (as that 
of Eq. ®). 

About the developped formalism, we are aware that the crucial point might be 
the practical realization of the generic linear quadratures such as in Eq. ([[]). Then, 
let us consider a practical implementation, in the optical domain. The quadrature 
of Eq. (|l|) could be experimentally accessible by using for example the squeezing 
pre-amplification (pre-attenuation) of a field mode which is going to be measured 
(a similar method in different context was discussed in Ref. [|34|])- In fact, let a 
be the signal field mode to be detected, when it passes through a squeezer it be- 
comes a s = a cosh s — a)e td sinh s, where s and 9 characterize the complex squeezing 



parameter ( = se p5| . Then, if we subsequently detect the field by using the bal- 
anced homodyne scheme, we get an output signal proportiopnal to the average of the 
following quadrature 

E^) = ^=(a s e~^ + a\e^) , (62) 

where is the local oscillator phase. When this phase is locked to that of the squeezer, 
such that = 9/2, Eq. (|62"D becomes 



E(4>) = (ae~ ie/2 [cosh s - sinh s] + a f e ie/2 [cosh s - sinh s]) , (63) 

which, essentially, coincides with Eq. ([!]), if one recognizes the independent parame- 
ters 

H = [cosh s — sinh s] cos(6'/2); v = [cosh s — sinh s] sin(0/2) . (64) 



The shift parameter 5 has not a real physical meaning, since it causes only a dis- 
placement of the distribution along the X line without changing its shape, as can be 
evicted from Eqs. (|4]) and (|5|). So, in a practical situation it can be omitted. To be 
more precise, the shift parameter does not play a real physical role in the measure- 
ment process, it has been introduced for formal completeness and it expresses the 
possibility to achieve the desired marginal distribution by performing the measure- 
ments in an ensemble of frames which are each other shifted; (related method was 
early discussed in Ref. [[36]]). In an electro-optical system this only means to have 
the freedom of using different photocurrent scales in which the zero is shifted by a 
known amount. 

7 Connection with measurements in homodyne to- 
mography 

At this point, a comparision with the usual tomographic technique, used in the ex- 
periments of the type of Ref. [17j], is useful. To this end we recall that in this case the 
timelike evolution of the system is brought about by the parameters changing, thus 
no explicit time dependence of w is needed. Furthermore, we note that a relation 
betweeen the density operator and the marginal distribution analogous to Eq. (|8|) 
can be derived starting from another operator identity such as 

J IX 

which, by the change of variables /i = -\/23«, v = \/23?a, becomes 

j 3= J- / dfidu Tr{pe-^}e^ = / dfxdu Tr{p e - t£ }e i£ . (66) 

The trace can be now evaluated using the complete set of eigenvectors {\x)} for the 
operator x, obtaining 

Tr{pe~"} = J dx w(x, fi, v)e~ ix (67) 
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then, putting this one into Eq. (|66|), we have a relation of the same form of Eq. (|8|) 
with the kernel given by 

K^ v = ^e' lx e ix = -Le^V^"^, (68) 

which is the same of Eq. (|S|) setting z — 1. It means that we now have only one 
particular Fourier component due to the particular change of variables (the most 
general should be z/j, = — i/23a and zv = a/23? a). 

In order to reconstruct the usual tomographic formula for the homodyne detection 



19| we need to pass in polar variables, i.e. fx = — rcos0, v = — rsin0, then 

x — > — rXcf, = — r[q cos 4> +psxn<p]. (69) 
Furthermore, indicating with x$ the eigenvalues of the quadrature x$, we have 

Tr{p e - l& } = Tr{pe 4ri >} = / dx^ w{x^ 4>)e irx « (70) 



p= d^dxtj, w(x^, (j))K^ (71) 



and thus, from Eq. (p 



with 



ks = — I dr re ir{x *-^ ] (72) 
2n J 

which is the same of Ref . [|B| . Substantially, the kernel of Eq. ( [72]) is given by the 
radial integral of the kernel of Eq. (|68D, and this is due to the fact that we go from 
a general transformation, with two free parameters, to a particular transformation 
(homodyne rotation) with only one free parameter, and then we need to integrate 
over the other one. This derivation follows Ref. p0| . 



8 Generating function for momenta 

Since the marginal distribution w(x, fi, v) has all the properties of the classical prob- 
ability density, one could calculate highest momenta for the shifted and sqeezed 
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quadrature x. We have by definition 

(x n ) = J x n w(x 1 /j,,iy)dx 1 (73) 

thus for the mean value (n = 1) 

(x) = J xw(x, /i, v)dx, (74) 

and for the quadrature dispersion one has 

&x = (x 2 ) — (x) 2 = J x 2 w(x, /i, v)dx — J xw(x, fj,,v)dx . (75) 

As in the standard probability theory [[37]], to calculate highest momenta for the 
shifted and squeezed quadrature one could introduces the generating function 

G(iX) = J2^-(x"). (76) 

Then the highest momenta are the coefficients of the Taylor series for the decompo- 
sition of the generating function with respect to the parameter (iX). We will express 
this generating function in terms of the Wigner function for the quantum system. 
Inserting Eq. (|5|) into Eq. ( [73] ) we have 

dqdpdkdx 

and inserting this one into Eq. © we arrive at 

G(iX) = I e -ik(*-w-»P)+iX* W ( qj p ) dqdpdkdx 
J (2tt) 2 

Now integrating, first over the quadrature variable x and then over the variable k, 

we get 

G(i\) = [ ^W(q,p)e iX ^ +up \ (79) 

This expression shows that the generating function for the quadrature highest mo- 
menta is determined by the Fourier components of the system Wigner function 

W F (a, b) = JL J W{q,p)e^ a+ipb dqdp (80) 
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(x n ) = / x n e- ik ^-^W(q, p) aq ^ ax (77) 



i. e. 



G(iX) = — W F (Xfi,Xu). (81) 

Thus, having the Wigner function of the quantum state and calculating its Fourier 
component, we may determine the generating function which depends on the extra 
parameters fi, v. On the other hand, since from Eq. fl30| ) we have the inverse Fourier 
transform 

W(q,p) = J W F (a,b)e-* qa - ipb dadb (82) 

we can express the Wigner function trough the generating function as 

W(q, P ) = e-^-^G(i)d^ (83) 

where we have taken A = 1 and integrated over the parameters fi and v on which the 
generating function depends. 

Hence we conclude that the quantum information about the state is completely 
contained in the expression for the generating function. It reflects the fact that 
measuring the shifted, rotated and squeezed quadrature we measure the momenta of 
the marginal distribution w(x,fi, u), and in fact we could reconstruct the generating 
function as a function of the extra parameters fi, v. Thus, the Wigner function of the 
system is obtained from Eq. (|83|). 

9 Conditional Probability 

The direct extension of classical probability concepts leads also to the conditional 
probability notion. Using the convention that x means the vector given by the quadra- 
ture variable x and the parameters ft and v, the joint probability u>(xx, ii; X2, £2) is 
defined as the probability to have X\ as result of the quadrature measurement at time 
t\ in the frame {fix, Ux} and X2 as result of the quadrature measurement at time £2 in 
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the frame {// 2 , ^2}- Then the conditional probability follows as 

u>(xi,£i;x 2 ,t 2 ) 



ltf(Xi,ti x 2 ,t 2 ) 



w(x 2 ,t 2 



M) 



As a consequence the Chapman-Kolmogorov equation [38] will be satisfied, i.e 



w(xi,tx) = / rf 3 x 2 u;(xi,ti|x 2 ,t 2 )w(x 2 ,t 2 ) 



55) 



so that the defined conditional probability can be interpreted as the propagator for 



the marginal distribution. The physical meaning of the real positive propagator (84]) 
is the following: it is the transition probability to go from the position x 2 in which 
the particle is situated at initial time t 2 i n the reference frame labeled by scaling and 
rotation parameters {fi 2 ,u 2 }, into the position x\ at the moment t\ in the reference 
frame labeled by the parameters {fii, 

We would remark that, even thought the stochastic process on which the marginal 
distribution depends is only one, we need to integrate also on the variables represent- 
ing the parameters since the same process may " come" from different frames. In fact 
really the normalization condition, as consequence of Eqs. @ and fl85l) , can be read 
as 

J dx 1 dx 2 w{x 1 ,t 1 \x2,t 2 )w(-x 2 ,t 2 ) = 1 . (86) 

In order to see the equation at which the conditional probability ( [84]) obeys, we 
insert Eq. (R3) into Eq. (BO), obtaining 



w(xi,ti|x 2 ,t 2 )(<9 tl + U(p,q)) = 5 (xi - x 2 )5(ti - h) 



17) 



that is the analogous of the differential Chapman-Kolmogorov equation [38 . 

As an application, let us consider the case of the driven harmonic oscillator for 
which, from Eq. (|3~4]) , we have 











ov\ a/ii 



9 \ 

+ f^g^rj ^(Xi,ti|x 2 ,t 2 ) = 5 3 (xi - X 2 )5(t! - h) 



2.3 



whose solution, for ti > t 2 , will be 



w(xi,ti|x 2 ,t 2 ) = 5(v 2 - /J,ism.[uj(t 1 -t2)]—ojv 1 cos[uj(t 1 — t 2 )]) 



x 5 (/ii cos[cj(ti — t 2 )} — <jJV\ cos[u;(ti — t 2 )] ~ ^2) 



/ 



/ . 



x 8 xi - x 2 — fii— ^{1 - cos[cu(ti - t 2 )]} — ^1— sin[u;(ti - t 2 )\ . 



10 



10 



19) 



Now, by means of Eqs. (85) and (|89| ) we may derive the solution of Eq. ( 34]) 
starting for example from an initial coherent condition characterized by go and po, 
i.e. Eq. ( pE2|) at t = 0, obtaining 

w(x, //, v) 



x exp 



^(/i 2 + uj 2 v 2 ) 

f f 
x — fi— -(1 — cos(a>£)) — sin(u;t) — go(yU cos(u;t) — ljv sm(ut)) 

CO 2 UJ 



(/^ + u;V) 



(90) 



where we have taken Xi = x, t\ = t and t 2 = 0. Of course, if we set / = and u = 1 
in Eq. (|9~0| ) we have again the solution fl42|). 

Finally, as special case of the propagator formula fl85|), we can consider the time 
evolution of the marginal distribution of Ref. pi 



w(xi, fix = cos 4>, v\ = sin (ft, t\] 



J d 3 x 2 w(xi,fii = cos0, V\ = sin0, ti|x 2 ,t 2 )w(x 2 ,t 2 ) 



(91) 



This could be useful as a connection between our formalism and the homodyne to- 
mography at different times. 
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10 Conclusions 



We have shown that it is possible to bring the quantum dynamics back to the classical 
description in terms of a probability distribution containing (over) complete informa- 
tion. The time evolution of a measurable probability for the discussed observables 
could be useful both for the prediction of the experimental outcomes at a given time 
and, as mentioned above, to achieve the quantum state of the system at any time. 
Furthermore, the symplectic transformation of Eq. ([[]) could be represented as a 
composition of shift, rotation and squeezing. So, we would emphasize that our pro- 
cedure allows to transform the problem of quantum measurements (at least for some 
observables) into a problem of classical measurements with an ensemble of shifted, 
rotated and scaled reference frames in the (classical) phase space. 

Quite generally physics distinguishes between the dynamical law and the state of 
a system. The state contains the complete statistical information about an ensemble 
of physical objects at a particular moment, while the dynamical law determines the 
change of the status quo at the next istant of time. But can we use the dynamical law 
to infer the state (for example) of a moving particle after position measurements have 
been performed? For istance, in molecular emission tomography |39| the quantum 
state of a mulecular vibration has been determined from its elongation encoded in the 
time-evolved fluorescence spectrum, while the usual standard tomography schemes 
ri7| have been restricted to harmonic oscillators or free particles for which one has 
a simple shearing or rotation in the phase space; however the developed formalism 
is able to infer the state of a particle moving in an arbitrary potential provided 



to have positionlike measurements in different frames (an analogous problem using 
nontomographic approach has been studied in Ref. f4]J). Of course, in some sit- 
uations the measurements of instantaneous values of the marginal distribution for 
different values of the parameters could be replaced by measuring the distribution 
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for these parameters which change in time. Such measurements may be consistent 
with the system evolution if the time variation of parameters is much faster than the 
natural evolution of the system itself. In this case the state of the system does not 
change during the measurement process and one obtains the instantaneous value of 
the marginal distribution and that of the Wigner function. 

Finally, we believe that our "classical" approach could be a powerful tool to in- 
vestigate complex quantum systems as for example chaotic systems in which the 
quantum chaos could be considered in a frame of equations for a real and positive 
distribution function. On the other hand, since the symplectic transformations are 
usually involved in the theory of special relativity, we could think to apply the devel- 
opped formalism for a relativistic formulation of the quantum measurement theory. 
These will be the subjects of future papers. 
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